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Abstract. The new generation of 3D hydrodynamical model atmospheres have 
been employed to study the impact of a realistic treatment of stellar convec- 
tion on element abundance determinations of globular cluster stars for a range of 
atomic and molecular lines. Due to the vastly different temperature structures in 
the optically thin atmospheric layers in 3D metal-poor models compared with cor- 
responding hydrostatic ID models, some species can be suspected to be hampered 
by large systematic errors in existing analyses. In particular, ID analyses based 
on minority species and low excitation lines may overestimate the abundances by 
> 0.3 dex. Even more misleading may be the use of molecular lines for metal-poor 
globular clusters. However, the prominent observed abundance (anti-)correlations 
and cluster variations are largely immune to the choice of model atmospheres. 
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1. Introduction mospheres are one-dimensional (ID, ei- 
ther plane-parallel or spherical), time- 
Determining stellar element abundances independent, static and assumed to fulfull 
play a crucial role in the efforts to improve hydrostatic equilibrium. Energy transport 
our understanding of formation and evo- by convection is approximated by the rudi- 
lution of globular clusters. The term ob- mentary mixing-length theory while oth- 
served abundances is somewhat of a mis- erwise radiative equilibrium is enforced, 
nomer however, since the chemical com- Furthermore, local thermodynamic equilib- 
position can not be inferred directly from rium (LTE) is normally assumed both for 
an observed spectrum. The obtained stel- the construction of the model atmospheres 
lar abundances are therefore never more and in the spectrum synthesis. It should 
trustworthy than the models of the stellar come as no surprise that abundance anal- 
atmospheres and the line formation pro- yses performed along these lines may well 
cesses employed to analyse the observa- contain significant systematic errors due to 
tions. Traditionally, abundance analyses of the adopted simplifications and approxima- 
late-type stars rely on a number of as- tions. 
sumptions, several of which are known to 

be of quite questionable nature. In stan- Perhaps the most severe shortcoming 

dard analyses, the employed model at- in standard analyses is the treatment of 
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convection. For late-type stars, the surface 
convection zone reaches the stellar atmo- 
sphere, which thereby directly affects the 
emergent spectrum. The solar granulation 
is the observational manifestation of con- 
vection: concentrated, rapid, cold down- 
drafts in the midst of broad, slow, warm 
upflows. Qualitatively similar granulation 
properties are expected in other solar-type 
stars, as indeed confirmed by 3D numeri- 
cal simulations (e.g. Nordlund & Dravins 
1990; Asplund et al. 1999; Asplund & 
Garcia Perez 2001; Allcndc Prieto et al. 
2002) and indicated by observed spectral 
line asymmetries. The up- and downflows 
have radically different temperature struc- 
tures (Stein & Nordlund 1998), which can 
not be approximated by normal theoreti- 
cal ID hydrostatic model atmospheres with 
different effective temperatures T e g (Fig. 
1). Because of the photospheric inhomo- 
geneities and the highly non-linear and 
non-local nature of spectrum formation, it 
is clear that no single ID model can be ex- 
pected to properly describe all aspects of 
what is inherently a 3D phenomenon (e.g. 
Asplund et al. 2003b). 

Here I will describe recent progress in 
developing 3D hydrodynamical model at- 
mospheres of late-type stars and their ap- 
plications to stellar abundance analyses, in 
particular for elements relevant for globular 
cluster studies. 

2. 3D hydrodynamical model 
atmospheres 

The 3D model atmospheres which form 
the basis of the abundance analyses pre- 
sented here have been computed with a 
3D, time-dependent, compressible, explicit, 
radiative-hydrodynamics code developed to 
study solar and stellar surface convection 
(Stein & Nordlund 1998). The hydrody- 
namical equations for conservation of mass, 
momentum and energy are solved on a 
Eulerian mesh with gridsizes of « 100 3 with 
explicit time-integration. The physical di- 
mensions of the grids are sufficiently large 
to cover many (> 10) granules simultane- 



ously in the horizontal direction and about 
13 pressure scale-heights in the vertical. In 
terms of continuum optical depth the simu- 
lations extend at least up to logTR, ss ~ —5 
which for most purposes are sufficient to 
avoid numerical artifacts of the open upper 
boundary on spectral line formation. The 
lower boundary is located at large depths to 
ensure that the inflowing gas is isentropic 
and featureless, while periodic horizontal 
boundary conditions are employed. 

In order to obtain a realistic atmo- 
spheric structure, it is crucial to have the 
best possible input physics, and properly 
account for the energy exchange between 
the radiation field and the gas. The adopted 
cquation-of-state is that of Mihalas et al. 
(1988), which includes the effects of ioniza- 
tion, excitation and dissociation. The con- 
tinuous opacities come from the Uppsala 
package (Gustafsson et al. 1975 and sub- 
sequent updates) while the line opacities 
are from Kurucz (1998, private communi- 
cation). The 3D radiative transfer is solved 
at each time-step under the assumptions 
of local thermodynamic equilibrium (LTE, 
S\ = B\) and opacity binning (Nordlund 
1982). The assignment of the original 2748 
wavelength points into the different opacity 
bins follows from detailed monochromatic 
radiative transfer calculations of the ID av- 
eraged atmospheric structure. The opacity 
binning thus includes the effects of line- 
blanketing in a manner reminiscent of opac- 
ity distribution functions. 

It is important to realise that the sim- 
ulations contain no free parameters which 
are tuned to improve the agreement with 
observations. The adoption of the numer- 
ical and physical dimensions of the simu- 
lation box is determined by practical com- 
putional time considerations, the need to 
resolve the most important spatial scales 
and the wish to place the artificial bound- 
aries as far as possible from the region of 
interest. It has been verified that the result- 
ing atmospheric structures are insensitive 
to the adopted effective viscosity at the cur- 
rent highest affordable numerical resolution 
(Asplund et al. 2000a). The input parame- 
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Fig. 1. The resulting temperature distribution in the upper part of 3D hydrodynamical 
convection simulations of the Sun (upper panel) and a metal-poor ([Fe/H] = —3) Sun 
(lower panel). Also shown are the predictions from the corresponding theoretical ID 
MARCS model atmospheres (dashed lines). The zero-point for the depth-scale corresponds 
roughly to the continuum optical depth unity. Note the much cooler 3D temperatures 
in the optically thin layers in the metal-poor star but the rough agreement at solar 
metallicity (see text for discussion). 
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ters discriminating different models are the 
surface gravity log g, metallicity [Fe/H] and 
the entropy of the inflowing material at the 
bottom boundary. The effective tempera- 
ture of the simulation is therefore a prop- 
erty which depends on the entropy struc- 
ture and evolves with time around its mean 
value following changes in the granulation 
pattern. 

The 3D hydrodynamical model atmo- 
spheres described above have sofar been 
performed for solar-type main sequence 
and subgiant stars (e.g. Nordlund & 
Dravins 1990; Asplund et al. 1999; Asplund 
& Garcia Perez 2001; Allende Prieto et al. 
2002). The most profound differences with 
predictions from theoretical ID hydrostatic 
models occur for low-metallicity stars, as 
first shown by Asplund et al. (1999). In ID. 
the presence of spectral lines causes sur- 
face cooling and backwarming (e.g. Mihalas 
1978), which translates to a shallow tem- 
perature gradient for metal-poor stars. In 
reality however, the radiative equilibrium 
which is enforced in ID model atmospheres 
is not necessary fulfilled. Instead the tem- 
perature in the optically thin atmospheric 
layers is determined by a competition be- 
tween two opposing effects: radiative heat- 
ing by absorption of photons in spectral 
lines and expansion cooling of overshoot- 
ing upflowing material due to the density 
stratification. At solar metallicity the two 
effects nearly balance, leaving the temper- 
ature close to the radiative equilibrium ex- 
pectations. At low metallicities, however, 
the lack of lines produces much less radia- 
tive heating and consequently the resulting 
temperature is much below the radiative 
equilibrium value (Asplund et al. 1999), as 
shown in Fig. □ At [Fe/H] = -3, the dif- 
ference between ID and 3D predictions can 
exceed 1000 K, which obviously can have a 
dramatic impact on spectral features sensi- 
tive to those cool layers, such as molecular 
and low excitation lines as well as minority 
species. 

Further details of the 3D hydrodynam- 
ical model atmospheres are available in 
Stein & Nordlund (1998), Asplund et al. 



(1999, 2000a,b) and Asplund & Garcia 
Perez (2001). 

3. 3D spectral line formation 

The 3D hydrodynamical model atmo- 
spheres described in the previous section 
form the basis for the 3D spectral line for- 
mation calculations. The original 3D hy- 
drodynamical simulations which includes 
deep, extremely optically thick layers are 
interpolated to a finer depth scale which 
extends only down to logTR OSS w 2.5 prior 
to the spectral line calculations for im- 
proved numerical accuracy. At the same 
time, the simulations are interpolated to 
a coarser horizontal grid to ease the com- 
putational burden. For abundance analy- 
sis purposes this procedure introduces non- 
noticable differences (< 0.01 dex for indi- 
vidual snapshots and less for temporal av- 
erages) . The concepts of micro- and macro- 
turbulence, which are introduced in ID 
analyses in order to account for the miss- 
ing line broadening, are not necessary in 
3D calculations with a self-consistent ac- 
counting of Doppler shifts arising from the 
convective motions. 

For the spectral line formation calcu- 
lations presented here the simplifying as- 
sumption of LTE has been made, implying 
that the level populations are determined 
by the Saha and Boltzmann distributions. 
With the source function thus known it is 
then straightforward, albeit computation- 
ally intensive, to solve the 3D radiative 
transfer equation for a number of simu- 
lation snapshots before spatial and tem- 
poral averaging of the resulting flux pro- 
files. All in all, a single temporally and 
spatially averaged flux profile in 3D cor- 
respond most of the time to N t * N x * 
N y * Aggies * N\ 10 8 ID radiative trans- 
fer calculations. In addition, each 3D pro- 
file is normally computed for at least three 
different abundances to enable interpola- 
tion to the requested line strength. Even 
then, such 3D LTE line calculations are 
achievable on current workstations thanks 
to efficient numerical algorithms. Recently, 
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Fig. 2. Schematic representation of the expected 3D LTE effects for different elements 
for T e ff = 5800 — 6200 K stars with [Fe/H] ~ —2; although not quantitatively investi- 
gated here stars with different parameters are expected to show qualitatively similar 3D 
effects. While many elements apparently are insensitive to the employment of 3D model 
atomspheres, some species can be quite vulnerable due to their temperature sensitivity, 
in particular molecular lines (e.g. CH, NH, OH) and low excitation lines. For example, 
while standard ID analyses yield accurate results for Fen lines, Fel lines tend to signifi- 
cantly overestimate the Fe abundances. Note that these estimates are based on 3D LTE 
calculations, while several species may well show significant departures from LTE. 



methods enabling even detailed 3D non- 
LTE line formation for large model atoms 
have been designed (Botnen & Carlsson 
1999) and applied to abundance analyses 
for Li (Kiselman 1997: Uitenbroek 1998; 
Asplund et al. 2003a) and O (Asplund et 
al. 2003b). The near future will doubtless 
see many more such studies for more ele- 
ments and 3D model atmospheres. 

In order to investigate the possible im- 
pact of the new generation of 3D hydro- 
dynamical model atmospheres on globular 
cluster research, I have selected a num- 
ber of species and lines which are often 
used to infer the chemical compositions of 
globular cluster stars: Li I, Ci, Oi, [Oi], 
Nai, Mgi, All, Si, Cai, Fei, Fen, Zni, 
Srn, Ban, CH, NH, and OH lines. The 



3D LTE line formation calculations were 
performed for solar-like (T g ~ 5800 K, 
loggf = 4.4) and turn-off stars (T g ~ 
6200 K, logg = 4.0) of varying metal- 
licities ([Fe/H] = 0.0,-1.0,-2.0,-3.0), as 
well as for a few 3D models correspond- 
ing to specific stars (Procyon, HD 140283, 
HD 84937, G64-12). Unfortunately, no 3D 
models are yet available for giants although 
work towards achieving this goal is on- 
going. As mentioned above, no microtur- 
bulence enters these 3D calculations. To 
facilitate an estimation of the 3D effects 
on derived abundances, corresponding cal- 
culations were carried out with the same 
code using ID MARCS model atmospheres 
(Asplund et al. 1997), adopting in all cases 
£turb = l.Okms -1 . 
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4. Results and discussion 

Fortunately, several elements appear little, 
if at all, affected by the employment of 
3D hydrodynamical model atmospheres in- 
stead of classical ID hydrostatic models, as 
shown in Fig. In particular, lines of Lil 
(in 3D non-LTE), Ci, Oi, Nai, Si, Fen, 
and Zni, as well as several Mgi and All 
lines, are essentially immune to the choice 
of model atmospheres (Aloge < 0.1 dex). 

The bad news, however, is that as ex- 
pected in metal-poor star minority species 
(e.g. Fei), low excitation (e.g. Mgi 517nm, 
All 394 nm, Srii 421 nm, Ban 455 nm) 
and all molecular lines tend to significantly 
over-estimate the abundances by > 0.3 dex 
when relying on ID model atmospheres. 
The reason for these differences is the much 
lower temperatures in the optically thin 
layers at low metallicities. It should be 
noted, however, that the use of 3D mod- 
els have recently caused a significant down- 
ward revision even of the solar C, N and O 
abundances (e.g. Asplund et al. 2003b) due 
mainly to the presence of temperature in- 
homogeneities in 3D. 

As a consequence of the above- 
mentioned findings, it is clear that globu- 
lar cluster metallicities should be based on 
Fen lines (see also Kraft & Ivans 2003). 
In addition, the observed O-Na, Al-Mg, 
C-N anti-correlations are robust features 
of globular clusters as are the large clus- 
ter abundance variations. While absolute 
abundances based on molecular lines can 
be highly misleading, isotopic abundance 
ratios are not significantly affected. 

It should be emphasized, however, that 
the above estimates are based only on 
3D LTE calculations. In some cases there 
may also be significant 3D non-LTE ef- 
fects, which may diminish the overall net 
effects of the new 3D model atmospheres 
for some species (e.g. Lil, Asplund et al. 
2003a) while aggrevate them for others 
(e.g. Oi, Asplund et al. 2003b). Detailed 
3D non-LTE calculations for additional ele- 
ments can be expected for the future, which 



would finally place the abundance analyses 
of late-type stars on a firm footing. 
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